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Abstract — In this paper, we study the problem of determining 
the size of battery storage used in grid-connected photovoltaic 
(PV) systems. In our setting, electricity is generated from PV 
and is used to supply the demand from loads. Excess electricity 
generated from the PV can be stored in a battery to be used later 
on, and electricity must be purchased from the electric grid if the 
PV generation and battery discharging cannot meet the demand. 
Due to the time-of-use electricity pricing, electricity can also be 
purchased from the grid when the price is low, and be sold back 
to the grid when the price is high. The objective is to minimize 
the cost associated with purchasing from (or selling back to) the 
electric grid and the battery capacity loss while at the same time 
satisfying the load and reducing the peak electricity purchase 
from the grid. Essentially, the objective function depends on the 
chosen battery size. We want to find a unique critical value 
(denoted as C^a) of the battery size such that the total cost 
remains the same if the battery size is larger than or equal to 
Cre(, and the cost is strictly larger if the battery size is smaller 
than Cref. We obtain a criterion for evaluating the economic 
value of batteries compared to purchasing electricity from the 
grid, propose lower and upper bounds on C^ef, and introduce 
an efficient algorithm for calculating its value; these results are 
validated via simulations. 

I. Introduction 

The need to reduce greenhouse gas emissions due to fossil 
fuels and the liberalization of the electricity market have led 
to large scale development of renewable energy generators 
in electric grids |r|. Among renewable energy technologies 
such as hydroelectric, photovoltaic (PV), wind, geothermal, 
biomass, and tidal systems, grid-connected solar PV continued 
to be the fastest growing power generation technology, with 
a 70% increase in existing capacity to 13GW in 2008 [2J. 
However, solar energy generation tends to be variable due to 
the diurnal cycle of the solar geometry and clouds. Storage 
devices (such as batteries, ultracapacitors, compressed air, and 
pumped hydro storage |[3]) can be used to i) smooth out the 
fluctuation of the PV output fed into electric grids ("capacity 
firming") ||2l, |]4|, ii) discharge and augment the PV output 
during times of peak energy usage ("peak shaving") ||5l, or 
iii) store energy for nighttime use, for example in zero-energy 
buildings and residential homes. 

Depending on the specific application (whether it is off- 
grid or grid-connected), battery storage size is determined 
based on the battery specifications such as the battery storage 
capacity (and the minimum battery charging/discharging time). 
For off-grid applications, batteries have to fulfill the following 
requirements: (i) the discharging rate has to be larger than 
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or equal to the peak load capacity; (ii) the battery storage 
capacity has to be large enough to supply the largest night time 
energy use and to be able to supply energy during the longest 
cloudy period (autonomy). The IEEE standard ||6| provides 
sizing recommendations for lead-acid batteries in stand-alone 
PV systems. In [7|, the solar panel size and the battery size 
have been selected via simulations to optimize the operation of 
a stand-alone PV system, which considers reliability measures 
in terms of loss of load hours, the energy loss and the 
total cost. In contrast, if the PV system is grid-connected, 
autonomy is a secondary goal; instead, batteries can reduce 
the fluctuation of PV output or provide economic benefits such 
as demand charge reduction, and arbitrage. The work in IH 
analyzes the relation between available battery capacity and 
output smoothing, and estimates the required battery capacity 
using simulations. In addition, the battery sizing problem has 
been studied for wind power applications f9l-fTTl and hybrid 
wind/solar power applications [12|-|14|. In [9J, design of a 
battery energy storage system is examined for the purpose 
of attenuating the effects of unsteady input power from wind 
farms, and solution to the problem via a computational proce- 
dure results in the determination of the battery energy storage 
system's capacity. Similarly, in [11 1, based on the statistics of 
long-term wind speed data captured at the farm, a dispatch 
strategy is proposed which allows the battery capacity to be 
determined so as to maximize a defined service lifetime/unit 
cost index of the energy storage system; then a numerical 
approach is used due to the lack of an expUcit mathematical 
expression to describe the lifetime as a function of the battery 
capacity. In |10|, sizing and control methodologies for a 
zinc-bromine flow battery-based energy storage system are 
proposed to minimize the cost of the energy storage system. 
However, the sizing of the battery is significantly impacted by 
specific control strategies. In llT2l . a methodology for calcu- 
lating the optimum size of a battery bank and the PV array 
for a stand-alone hybrid wind/PV system is developed, and a 
simulation model is used to examine different combinations 
of the number of PV modules and the number of batteries. 
In ||T3l . an approach is proposed to help designers determine 
the optimal design of a hybrid wind-solar power system; 
the proposed analysis employs linear programming techniques 
to minimize the average production cost of electricity while 
meeting the load requirements in a reliable manner. In |14|, 
genetic algorithms are used to optimally size the hybrid system 
components, i.e., to select the optimal wind turbine and PV 
rated power, battery energy storage system nominal capacity 
and inverter rating. The primary optimization objective is the 
minimization of the levelized energy cost of the island system 
over the entire lifetime of the project. 
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In this paper, we study the problem of determining the 
battery size for grid-connected PV systems. The targeted 
applications are primarily electricity customers with PV arrays 
"behind the meter," such as residential and commercial build- 
ings with rooftop PVs. In such applications, the objective is to 
reduce the energy cost and the loss of investment on storage 
devices due to aging effects while satisfying the loads and 
reducing peak electricity purchase from the grid (instead of 
smoothing out the fluctuation of the PV output fed into electric 
grids). Our setting is shown in Fig.[T] Electricity is generated 
from PV panels, and is used to supply different types of 
loads. Battery storage is used to either store excess electricity 
generated from PV systems for later use when PV generation 
is insufficient to serve the load, or purchase electricity from the 
grid when the time-of-use pricing is lower and sell back to the 
grid when the time-of-use pricing is higher. Without a battery, 
if the load was too large to be supplied by PV generated 
electricity, electricity would have to be purchased from the 
grid to meet the demand. Naturally, given the high cost of 
battery storage, the size of the battery storage should be chosen 
such that the cost of electricity purchase from the grid and the 
loss of investment on batteries are minimized. Intuitively, if 
the battery is too large, the electricity purchase cost could be 
the same as the case with a relatively smaller battery. In this 
paper, we show that there is a unique critical value (denoted 
as Cjgf, refer to Problem [T]) of the battery capacity such that 
the cost of electricity purchase and the loss of investments on 
batteries remains the same if the battery size is larger than or 
equal to C^^j and the cost is strictly larger if the battery size 
is smaller than C^^f. We obtain a criterion for evaluating the 
economic value of batteries compared to purchasing electricity 
from the grid, propose lower and upper bounds on C^^^ given 
the PV generation, loads, and the time period for minimizing 
the costs, and introduce an efficient algorithm for calculating 
the critical battery capacity based on the bounds; these results 
are validated via simulations. 

The contributions of this work are the following: i) to the 
best of our knowledge, this is the first attempt on determining 
the battery size for grid-connected PV systems based on a 
theoretical analysis on the lower and upper bounds of the 
battery size; in contrast, most previous work are based on 
simulations, e.g., the work in f8l- lfT2l ; ii) a criterion for eval- 
uating the economic value of batteries compared to purchasing 
electricity from the grid is derived (refer to Proposition |4] and 
Assumption |2j, which can be easily calculated and could be 
potentially used for choosing appropriate battery technologies 
for practical applications; and iii) lower and upper bounds on 
the battery size are proposed, and an efficient algorithm is 
introduced to calculate its value for the given PV generation 
and dynamic loads; these results are then validated using sim- 
ulations. Simulation results illustrate the benefits of employing 
batteries in grid-connected PV systems via peak shaving and 
cost reductions compared with the case without batteries (this 
is discussed in Section IVlB). 

'Note that solar panels and batteries both operate on DC, while the grid and 
loads operate on AC. Therefore, DC-to-AC and AC-to-DC power conversion 
is necessary when connecting solar panels and batteries with the grid and 
loads. 



The paper is organized as follows. In the next section, we lay 
out our setting, and formulate the storage size determination 
problem. Lower and upper bounds on C^^f are proposed 
in Section |III1 Algorithms are introduced in Section |IV] to 
calculate the value of the critical battery capacity. In SectionlVl 
we validate the results via simulations. Finally, conclusions 
and future directions are given in Section [Vl] 

II. Problem Formulation 

In this section, we formulate the problem of determining 
the storage size for a grid-connected PV system, as shown in 
Fig. [T] We first introduce different components in our setting. 

A. Photovoltaic Generation 

We use the following equation to calculate the electricity 
generated from solar panels: 

Ppv(t) = Gm{t) xSxr] , (1) 

where 

• GHI (Wm~^) is the global horizontal irradiation at the 
location of solar panels, 

> S (m^) is the total area of solar panels, and 

• 77 is the solar conversion efficiency of the PV cells. 
The PV generation model is a simplified version of the one 
used in (15] and does not account for PV panel temperature 
effects 

B. Electric Grid 

Electricity can be purchased from or sold back to the grid. 
For simplicity, we assume that the prices for sales and pur- 
chases at time t are identical and are denoted as Cg{t){$/Wh). 
Time-of-use pricing is used (Cg(<) > depends on t) because 
commercial buildings with PV systems that would consider 
a battery system usually pay time-of-use electricity rates. 
In addition, with increased deployment of smart meters and 
electric vehicles, some utility companies are moving towards 
residential time-of-use pricing as well; for example, SDG&E 

-Note that our analysis on determining battery capacity only relies on 
Ppv(t) instead of detailed models of the PV generation. Therefore, more 
complicated PV generation models can also be incorporated into the cost 
minimization problem as discussed in Section IlllF. 



Load 




Electric Grid Battery 



Fig. 1. Grid-connected PV system with battery storage and loads. 
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(San Diego Gas & Electric) has the peak, semipeak, offpeak 
prices for a day in the summer season [TE], as shown in Fig. [3] 
We use Pg{t){W) to denote the electric power exchanged 
with the grid with the interpretation that 

• Pg{t) > if electric power is purchased from the grid, 

• Pg{t) < if electric power is sold back to the grid. 

In this way, positive costs are associated with the electricity 
purchase from the grid, and negative costs with the electricity 
sold back to the grid. In this paper, peak shaving is enforced 
through the constraint that 

P,it) < D , 

where _D is a positive constant. 

C. Battery 
A battery has the following dynamic: 

dEB{t) 



dt 



Psit) , 



(2) 



where EB{t){Wh) is the amount of electricity stored in the 
battery at time t, and PB{t){W) is the charging/discharging 
rate; more specifically, 

• Psit) > if the battery is charging, 
> and Pb (t) < if the battery is discharging. 
For certain types of batteries, higher order models exist (e.g., 
a third order model is proposed in ifTTI . ifTSi ). 

To take into account battery aging, we use C{t){Wh) to de- 
note the usable battery capacity at time t. At the initial time to, 
the usable battery capacity is Ci-ef, i.e., C{to) ~ Cref > 0. The 
cumulative capacity loss at time t is denoted as AC{t){Wh), 
and AC (to) = 0. Therefore, 

C{t) = C,-ef - AC{t) . 

The battery aging satisfies the following dynamic equation 

dAC{t) _ i-ZPB{t) if PB{t) < 
10 



dt 



otherwise 



(3) 



where Z > is a constant depending on battery technologies. 
This aging model is derived from the aging model in |5| 
under certain reasonable assumptions; the detailed derivation 
is provided in the Appendix. Note that there is a capacity loss 
only when electricity is discharged from the battery. Therefore, 
AC(t) is a nonnegative and non-decreasing function of t. 
We consider the following constraints on the battery: 

i) At any time, the battery charge EB{t) should satisfy 

Q<EB{t)<C{t)^a,i~AC{t) , 

ii) The battery charging/discharging rate should satisfy 

-PBmin < Psit) < -Pfimax , 

where PBmin < 0, — PBmin IS the maximum battery 
discharging rate, and Pfimax > is the maximum battery 
charging rate. For simplicity, we assume that 

P -_P _ C{t) _ C-ef-ACft) 

where constant Tc > is the minimum time required to 
charge the battery from to C{t) or discharge the battery 
from C(t) to 0. 



D. Load 

P\o-dd{t){W) denotes the load at time t. We do not make 
explicit assumptions on the load considered in Section |lll] 
except that Pioad(i) is a (piecewise) continuous function. In 
residential home settings, loads could have a fixed schedule 
such as lights and TVs, or a relatively flexible schedule 
such as refrigerators and air conditioners. For example, air 
conditioners can be turned on and off with different schedules 
as long as the room temperature is within a comfortable range. 



E. Converters for PV and Battery 

Note that the PV, battery, grid, and loads are all connected 
to an AC bus. Since PV generation is operated on DC, a DC- 
to-AC converter is necessary, and its efficiency is assumed to 
be a constant rjpy satisfying 

< ?7pv < 1 . 

Since the battery is also operated on DC, an AC-to-DC 
converter is necessary when charging the battery, and a DC- 
to-AC converter is necessary when discharging as shown in 
Fig. [T] For simplicity, we assume that both converters have 
the same constant conversion efficiency tje satisfying 

< r?B < 1 . 

We define 



p mPBit) iiPB{t)<o 

PBC[t) = < p^(t) 

—s±J- otherwise , 

I. '?B ' 

In other words, PBc{t) is the power exchanged with the AC 
bus when the converters and the battery are treated as an entity. 
Similarly, we can derive 



( PbcW 

[VsPBcit) otherwise, 



if Pfic(t) < 



Note that rj^ is the round trip efficiency of the battery storage. 



E Cost Minimization 

With all the components introduced earlier, now we can 
formulate the following problem of minimizing the sum of 
the net power purchase cosj| and the cost associated with the 
battery capacity loss while guaranteeing that the demand from 



'Note that the net power purchase cost include the positive cost to purchase 
electricity from the grid, and the negative cost to sell electricity back to the 
grid. 
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loads and the peak shaving requirement are satisfied: 

i-to+T 

min / Cg(T)Pg(T)dT + i^AC(to+T) 

S.t. 77pv^pv(0 + Pg{t) = ftc(i) + Plo,6{t) , 



(6) 



dt 



= PB{t) 



dACjt) _ I -ZPeit) if Pij(t) < 

dt 1 otherwise , 

0<i;B(i) <C,ef-AC(i) , 
EB{to) = ,ACito) = , 

- AC{t) 



Pi 



Bmax — ^Bmin 



Ps(<) = < 

I 7?B-PBc(i) otherwise , 

Pg(t) < , 



(7) 



where to is the initial time, T is the time period considered 
for the cost minimization, K{$/Wh) > is the unit cost for 
the battery capacity loss. Note that: 

i) No cost is associated with PV generation. In other words, 
PV generated electricity is assumed free; 

ii) KAC {to + T) is the loss of the battery purchase invest- 
ment during the time period from to to tg + T due to the 
use of the battery to reduce the net power purchase cost; 

iii) Eq. (|6]l is the power balance requirement for any time 

[ta,to + T]; 

iv) the constraint Pg{t) < D captures the peak shaving 
requirement. 

Given a battery of initial capacity Cief, on the one hand, if 
the battery is rarely used, then the cost due to the capacity 
loss KAC {to + T) is low while the net power purchase 
cost J^"'^'^ Cg{T)Pg(T)dT is high; on the other hand, if the 
battery is used very often, then the net power purchase cost 
fto^ ^&i'^)P%i'^)'^'^ while the cost due to the capacity 

loss KAC{to + T) is high. Therefore, there is a tradeoff on 
the use of the battery, which is characterized by calculating 
an optimal control policy on Pb [t) , Pg [t) to the optimization 
problem in Eq. (|7]l- 

Remark 1 Besides the constraint Pg (t) < D, peak shaving is 
also accomplished indirectly through dynamic pricing. Time- 
of-use price margins and schedules are motivated by the 
peak load magnitude and timing. Minimizing the net power 
purchase cost results in battery discharge and reduction in grid 
purchase during peak times. If, however, the peak load for a 
customer falls into the off-peak time period, then the constraint 
on Pg{t) limits the amount of electricity that can be purchased. 
Peak shaving capabilities of the constraint Pg{t) < D and 
dynamic pricing will be illustrated in Section |V]B. ■ 

G. Storage Size Determination 
Based on Eq. (|6]l, we obtain 

Pg{t) = Pload(t) - ?7pvPpv(i) + Pucit) . 



Let u{t) — Pbc(0' then the optimization problem in Eq. (|7|l 
can be rewritten as 

min/;°+^ Cg(r)(Pi„,d(r) - 7]pvPpv(^) + <'r))dr + KAC{to + T) 



s.t. 



dt I r]Bu{t) otherwise 



u(t) 

otherwise , 



dACjt) _ J -Z^ if u{t) < 
di ^ 



EB{t)>0 ,EB{to) = ,AC{to)^0 , 

Eeit) + AC{t) < C,,f , 

r]Bu{t)T, + AC{t) < Crf if u{t) > 0, 

- + AC{t) < Cref if u{t) < 0, 

Pio.d{t) - V^,P^,{t) + u{t) < D . 



(8) 



Now it is clear that only u{t) is an independent variable. We 
define the set of feasible controls as controls that guarantee 
all the constraints in the optimization problem in Eq. dHJ. 
Let J denote the objective function 

min„ Cs{r){Pio-Ar) " V-Ppv(^) + ■4r))dT + KAC{to + T) . 

If we fix the parameters to,T, K, Z, Tc, and Z?, J is a function 
of Cief, which is denoted as J(Cref). If we increase Cref, 
intuitively J will decrease though may not strictly decrease 
(this is formally proved in Proposition |2]i because the battery 
can be utilized to decrease the cost by 

i) storing extra electricity generated from PV or purchasing 
electricity from the grid when the time-of-use pricing is 
low, and 

ii) supplying the load or selling back when the time-of-use 
pricing is high. 

Now we formulate the following storage size determination 
problem. 

Problem 1 (Storage Size Determination) Given the opti- 
mization problem in Eq. dHJ with fixed to , T, K, Z, Tc, and 
D, determine a critical value C^^f > such that 

. VGef < C^^f, J(aef) > J(Q'ef)' ^nd 
. VGef > C^^„ J(aef) = J(C;,f). 

One approach to calculate the critical value C^^f is that we 
first obtain an explicit expression for the function J(Cief) by 
solving the optimization problem in Eq. (O and then solve 
for C^^f based on the function J. However, the optimization 
problem in Eq. dHJ is difficult to solve due to the nonlinear 
constraints on u{t) and Esit), and the fact that it is hard to 
obtain analytical expressions for Pioad(i) and Ppv(i) in reality. 
Even though it might be possible to find the optimal control 
using the minimum principle |19|, it is still hard to get an 
explicit expression for the cost function J. Instead, in the next 
section, we identify conditions under which the storage size 
determination problem results in non-trivial solutions (namely, 
Cj^f is positive and finite), and then propose lower and upper 
bounds on the critical battery capacity C^^f. 
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III. Bounds on C^^^ 

Now we examine the cost minimization problem in Eq. ([8]). 
Since Pioad(i) — '7pv^pv(i) -^u{t) < D, or equivalently, u{t) < 
D + rij,yPpy{t) — Pioad(i)' is ^ constraint that has to be satisfied 
for any t e [to,to + T], there are scenarios in which either 
there is no feasible control or u{t) — for t E [to, to + T]. 
Given Ppv{t), Pioad{t), and D, we define 

51 = {te [to, to + T]\D + vPpv(t) - Pioad(i) < 0} , (9) 

52 = {te [to, to+T]\D + vPpv(i) - fload(i) = 0} , (10) 

53 = {te[to,to+T]\ D + vPpv(i) - fload(i) > 0} . (11) 

Note tha0 S"! © 5*2 ® 5*3 = [to,to + T]. Intuitively, Si is the set 
of time instants at which the battery can only be discharged, 
S2 is the set of time instants at which the battery can be 
discharged or is not used (i.e., u{t) — 0), and ^3 is the set of 
time instants at which the battery can be charged, discharged, 
or is not used. 

Proposition 1 Given the optimization problem in Eq. (O, if 

i) to e Si, or 

ii) 5*3 is empty, or 

iii) to ^ (or equivalently, 6 S'2 U S3), S3 is nonempty, 
5*1 is nonempty, and 3ti G Si, Vta S S3, ti < 13, 

then either there is no feasible control or u{t) — for t G 
[to,to+T]. 

Proof: Now we prove that, under these three cases, either 
there is no feasible control or u{t) = for t e [to,to + T]. 

i) If to £ Si, then u{to) <D + rj^^Pp^to) - Pioad(io) < 0. 
However, since Esito) = 0, the battery cannot be 
discharged at time to- Therefore, there is no feasible 
control. 

ii) If 5*3 is empty, it means that u{to) < D + »ypvPpv(^o) — 
-Pioad(^o) £ for any t € [to,to + T], which implies that 
the battery can never be charged. If Si is nonempty, then 
there exists some time instant when the battery has to be 
discharged. Since Esito) — and the battery can never 
be charged, there is no feasible control. If 5*1 is empty, 
then u{t) = Q for any t E [to,to+ T] is the only feasible 
control because £'5(^0) — 0. 

iii) In this case, the battery has to be discharged at time ti, 
but the charging can only happen at time instant is G S3. 
If Vt3 G 5*3, 3ti G 5*1 such that ti < ^3, then the battery 
is always discharged before possibly being charged. Since 
Esito) — 0, then there is no feasible control. 

■ 

Note that if the only feasible control is u{t) — for t G 
[to, to + T], then the battery is not used. Therefore, we impose 
the following assumption. 

Assumption 1 In the optimization problem in Eq. (HJ, to G 

5*2 U 53, S3 is nonempty, and either 

• 5*1 is empty, or 

• 5*1 is nonempty, but Vti G ^i, 3t3 G ^3, ts < ti, 
where Si,S2, S3 are defined in Eqs. (|9]), ([T0|, (fTTT) . 



Given Assumption [U there exists at least one feasible 
control. Now we examine how J(Cief) changes when Cret 
increases. 

Proposition 2 Consider the optimization problem in Eq. ^ 
with fixed tQ,T,K,Z,Tc, and D. If C^f < C^f, then 
J(Ci,.) > JiCl,). 

Proof: Given C^^f, suppose control u^{t) achieves the 
minimum cost J{C^^f) and the corresponding states for the 
battery charge and capacity loss are Eg{t) and AC^(t). Since 

Elit) + AC\t)<Ci,<Cl, , 
mu\t)Tc + AC^t) < Cif < Cli if u\t) > 0, 
u\t), 



-T, + ACHt) < < if u^t) < 0, 



PioAt) - Vp^Pp^it) + u\t) < D , 

u^{t) is also a feasible control for problem ^ with C^^j, and 
results in the cost J{C^^f). Since J{C'^^f) is the minimum cost 
over the set of all feasible controls which include u^{t), we 
must have J{C^^f) > J(Cjef)- ■ 
In other words, J is non-increasing with respect to the 
parameter Cief, i.e., J is monotonically decreasing (though 
may not be strictly monotonically decreasing). If Cref — 0, 
then < AC(t) < Cref = 0, which implies that u{t) = 0. In 
this case, J has the largest value 

Jmax = J(0)= /""^ Cg(T)(P,oad(T)-vPpv(T))dT . (12) 

Jto 

Proposition |2] also justifies the storage size determination 
problem. Note that the critical value C^^^ (as defined in 
Problem [T]) is unique as shown below. 

Proposition 3 Given the optimization problem in Eq. (HJ with 
fixed to, T, K, Z, Tc, and D, C^^f is unique. 

Proof: We prove it via contradiction. Suppose C^^f is not 
unique. In other words, there are two different critical values 
C^gf and C^gf . Without loss of generality, suppose C^^j < C^^j. 
By definition, J(C^ef) > ^i^r^) because C^^^ is a critical 
value, while J{C^^f) = J{C^^) because C^^j is a critical value. 
A contradiction. Therefore, we must have C^^j = C^^j. ■ 
Intuitively, if the unit cost for the battery capacity loss K is 
higher (compared with purchasing electricity from the grid), 
then it might be preferable that the battery is not used at all, 
which results in C^^f = 0, as shown below. 

Proposition 4 Consider the optimization problem in Eq. ([8j 
with fixed to , T, K, Z, Tc, and D under Assumption [T] 

i) If 



K > 



(maxt Cg{t) - mint Cg{t))riB 



1 B means C = AUB and AnB - 



then J(Cref) = Jmax, which implies that C^^j = 0; 
ii) if 

(maxt Cg(t) - mint Cg(t) )?7b 
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then 



KZ 

J(Ci-ef) >Jmax " (max Cg(t) - min Cg(t) )x 

T X (I? + max(77pvFpv(i) - Fioad(t))) , 



where maxf and mint are calculated for i e [tp , to + . 
Proof: The cost function can be rewritten as 



^(Cref) — 



Cg(T)u(T)dT + KAC{to + T) 



max ^ / 

to 

Jmax + / Cg(T)u(r)(ir+ 



Jt„ Vb 



u(t)<0 



dT 



where 



and 



to+T 



Cg(T)"('r)l«(r)>oC'T , 



J- = / (Cg(T) )u{T)\u{r)<odT . 

Jto VB 

Note that J+ > because the integrand Cg(T)M(T)|„(7.')>o is 
always nonnegative. 

i) We first consider K > ("^'^x, (t)-min, c,(t))^B ^ ^j^^^^ 
are two possibilities: 



. K > 



(maxt Cg(t))7)B 



, or equivalently, maXf Cg(t) < 

Thus, for any t G [to, to + T], Cg{t) - ^ < 0, which 
imphes that J_ > 0. Therefore, J(Cref) > Jmax- 

(maxt Cg{t) — mint Cg(t))r)B < ^ < (maxt Cg{t))7]B 

Let ^1 = maxt Co(t) - — , then Ai > 0. 
Let ^2 = mintCg(<), then ^2 > 0. Since 

(max,C,(t)-min,C,(t))„B < K , A2 > A^. NoW WC 
t +T 

have J+ > A2 /J M(T)|„(^)>odT, and J_ > 
Ai //^°^'^M(T)|„(^)<odT. Therefore, 

j-to+T 

JlCref) >Jn,ax + ^2 / "('t) | n(r) >0C?'^+ 

■Jto 

rto+T 

^1 / u(T)|„(T)<0C^'r 



rto+T 

= Jmax + {A2 ~ Ai) u(t)|„(^) 
/•to+T 

Ai / u(r)dT 

"'to 

(■to+T 

> Jmax + Ai u{T)dT . 

Jto 

Since Esito +T)= Eg (to) + j'"'^'^ PB{T)dT > and 



dT- 



Esito) — 0, we have 

/•to+T 
< / PB{T)dT 
Jto 

ta+T /•to+T 

PB{T)\ps(r)>odT + / PB{T)\p^(^r)<odT 

3 "'to 

to+T /•to+T ^(^^^ 

Jto VB 
rto+T /•to+T 

< / u{T)\u{r)>odT + / M(T)|„(r)<orfT 



to+T 



u{T)dT 



Therefore, we have J(Cref) > Jm 



(13) 



In summary, if K > ("''^x, c,(t)-min, c,(t))„B ^ j^^^^ 

J{Crei) > Jmax, which impHes that J(Cref) > J{0)- Since 
J(C'ref) is a non-increasing function of Cref, we also have 
J{Ccef) < J(0). Therefore, we must have J(Cref) = ^(0) = 
Jmax for Cref > 0, and Cj'gf = by definition. 



(maxt Cg(t)— mint Cg(t))i7B 



ii) We now consider K < ^ 
implies that Ai > A2 > 0. Then 



which 



rto + i 

J(Cref) >Jmax +^2 / u{t)\ 

m(t)>0 

Jto 

rto+T 

M / u{T)\u(r)<odT 



dT+ 



to+T 



^Jm^LX + A2 u{T)dT+ 
Jto 
rto+T 

{Ai ~ A2) / u(T)U(r)<0C^T • 

"'to 

Since J^"^'^ u(T)dT > as argued in the proof to i), 

J(aef) > Jmax + (^i - A2) //j^^ ?/(t) |„(,)<odT. NoW 

we try to lower bound J^"^"^ u(T)\u(^^^^odT. Note that 
u{T)dT > (as shown in Eq. ( fTSl )) implies that 

to+T /•to+T 

u{T)\u(T)<odr > - / u{T)\u{r)>adT . 

3 "'to 

Given Assumption [T] ^3 is nonempty, which implies that 

D + maxt(77pvPpv(t) - Pioad(t)) > 0. Since u{t) < D + 
f?pv-Ppv(i) - P\oRd{t) < D + maxt(?7pvPpv(0 - Pioad(i)), 
Jt°^^ uiT)\^(r)>odT < Tx (i^ + maxt(vPpv(t)-Pioad(t))), 
which implies that 

/,^°+^ «(T)|„(,)<odT > -T X (D + maxi(vPpv(t) - P|oad(t))) • (14) 

In summary, J(Cref) > Jmax ~ {Ai ~ A2) x T x {D + 
maxf (?7pvPpv(t) — ^'ioad(i))), which proves the result. ■ 

Remark 2 Note that i) in Proposition H] provides a criterion 
for evaluating the economic value of batteries compared to 
purchasing electricity from the grid. In other words, only if 
the unit cost of the battery capacity loss satisfies 



K < 



(maxt Cg(t) - mint Cg{t))T]B 



7 



it is desirable to use battery storages. This criterion depends on 
the pricing signal (especially the difference between the max- 
imum and the minimum of the pricing signal), the conversion 
efficiency of the battery converters, and the aging coefficient 
Z. ■ 

Given the result in Proposition H] we impose the following 
additional assumption on the unit cost of the battery capacity 
loss to guarantee that Cj^j is positive. 

Assumption 2 In the optimization problem in Eq. ([8]l, 

(maxt Cg(t) - mint Cg{t))vB 
Z 

In other words, the cost of the battery capacity loss during 
operations is less than the potential gain expressed as the 
margin between peak and off-peak prices modified by the 
conversion efficiency of the battery converters and the battery 
aging coefficient. Since J(Cief) is a non-increasing function of 
Cief and lower bounded by a finite value given Assumptions [T] 
and 121 the storage size determination problem is well defined. 
Now we show lower and upper bounds on the critical battery 
capacity in the following proposition. 

Proposition 5 Consider the optimization problem in Eq. ^ 
with fixed to, T, K, Z, Tc, and D under Assumptions [T] and |2] 
Then C;^^ < C^^f < C^^f, where 

= max(— X (max(Pioad(t) - ??pvPpv(t)) - D), 0) , 

and 



ZT 
'7b ' 



ryeT) x (L»+max(r/pvPpv(t)-/ioad(t)) 



Proof: We first show the lower bound via contradic- 
tion. Without loss of generality, we assume that x 
(maxt(Pioad(0 — '7pv^pv(i)) — D) > (because we require 
C,^f > 0). Suppose 

T 

ref ^ref 



^ X (max(Pioad(t) - 77pvPpv(0) - D) 



or equivalently. 



D < + max(Pioad(i) - r?pvPpv(t)) • 

Therefore, there exists ti g [to, to + T] such that D < 



--Pioad(ii)-7?pv-Ppv(ii)- Since u{ti) < £'-Pioad(ii)- 



Vp^PpAti) < < 0, we have 

„ u{ti) D + 77pvPpv(il) - fioadCii) 

^B[ti) — S 

CI, , aef-AC(ti) 



< 



< 



The implication is that the control does not satisfy the dis- 
charging constraint at <i. Therefore, C^^^ < C^^f. 

To show C^^f < C^*f, it is sufficient to show that if Cref > 
^ref' electricity that can be stored never exceeds C"Jj. 

Note that the battery charging is limited by Pfimax, i.e., 
-Pb(^) ~ T]Bu{t) < Psmax. Now wc try to lower bound Psmax. 

a-ef - AC{t) ^ Gef - AC(io + T) 



Pr 



> 



in which the second inequality holds because AC{t) is a 
non-decreasing function of t. By applying the aging model 
in Eq. ([3]l, we have 

Bmax ^ m — 



'-\-u,(T)<od.T 



Using Eq. (fT4] l. we have 

Cref - ^ X (D + maXt(?7pvPpv(0 - Pload(i))) 



Since Cref > C"gf, we have 



^ (max(i7Brc + |3:,77BT)-|^)x(D+maxt(»7p,Pp,(f)-Pioad(t))) 
-^Bmax _ 

>Vb{D + max(77pvPpv(t) - Pioad(i))) ■ 

Therefore, u{t) < D + ijpyPpy{t) - Pioad(t) < D + 
max((7ypvPpv(t) — Pioad(O) — "^^7^' which impHes that 
Psit) = rjBU^t) < Pfimax. Thus, the only constraint on u 
related to the battery charging is u{t) < I? + max((?7pvPpv(t) — 
Pioad(0)- In turn, during the time interval [to, to + T], the 
maximum amount of electricity that can be charged is 

fto+T nto+T 

PB{T)\ps(T)>odT ^ / ■m'^{T)[u(T)>odT 
'to Jto 

<77bT x{D + max(77pvPpv(t) - Pioad(O)) , 

which is less than or equal to CJ^Jj. In summary, if Cref > C^Jj, 
the amount of electricity that can be stored never exceeds C, 



ref 



and therefore, C^^^ < C^^f. ■ 

Remark 3 As discussed in Proposition [l] if 5*3 is empty, or 
equivalently, (i:) + maxt(?7pvPpv(t) -Pioad(t)) < 0, then C^^^ < 
0, which implies that C^^f = 0; this is consistent with the result 
in Proposition [T] when 5*1 is empty. ■ 

IV. Algorithms for Calculating C^^^ 

In this section, we study algorithms for calculating the 
critical battery capacity C^^j. 

Given the storage size determination problem, one approach 
to calculate the critical battery capacity is by calculating 
the function J(Cref) and then choosing the Cref such that 
the conditions in Problem [T] are satisfied. Though Cref is a 
continuous variable, we can only pick a finite number of Cref 's 
and then approximate the function J(Cref). One way to pick 
these values is that we choose Cef from C^Jj to C^'^j with the 
step siz^ — Tcap < 0. In other words. 



CI 



■ef ■ 



Cub 
ref 



cap ; 

y~nib ^Ib 

|-_rEf — [rf-| Suppose the picked 



wher^i = 0,l,...,LandL 
value is C^^f, and then we solve the optimization problem in 
Eq. (O with Cjgj. Since the battery dynamics and aging model 
are nonlinear functions of u{t), we introduce a binary indicator 
variable Iu(t)^ in which I^^t) = if u{t) > and Iu{t) = 1 if 

'Note that one way to implement the critical battery capacity in practice is 
to connect multiple identical batteries of fixed capacity Cflxcd in parallel. In 
this case, Tcap can be chosen to be Cflxcd- 

^The ceiling function \x~\ is the smallest integer which is larger than or 
equal to x. 
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Algorithm 1 Simple Algorithm for Calculating C[ 



ref 



Algorithm 2 Efficient Algorithm for Calculating C[ 



Input: The optimization problem in Eq. ^ with fixed 
to,T, K, Z,Tc, D under Assumptions [T] and |2l the calculated 
bounds C,'ef, C"ef, and parameters 5t, reap, Ttost 

Output: An approximation of C^ef 



1: Initialize C^f 

-c^.f-c'';,-,. 



■cf 



j X Trap, where i = 0, 1, L and L 



r- 



T"c:ap ' 

2: for i = 0, 1, L do 

3: Solve the optimization problem in Eq. ^ with C,'ef, and obtain 



J(c;e,); 

if i > 1 and J{Cl,f) - J(Cjif) > r,.osi then 
Set Cjef — C,'e}, and exit the for loop; 

end if 
end for 
Output C^^f. 



u{t) < 0. In addition, we use 5t as the sampling interval, and 
discretize Eqs. (|2j and ^ as 



AC{k- 



Eeik) + PB{k)St , 

jAC{k) - ZPB{k)St iiPB{k)<0 
|AC(fc) otherwise. 



With the indicator variable Iu{t} and the discretization of 
continuous dynamics, the optimization problem in Eq. ^ 
can be converted to a mixed integer programming problem 
with indicator constraints (such constraints are introduced in 
CPlex f20|), and can be solved using the CPlex solver [i20l 
to obtain J{Cl.^f). In the storage size determination problem, 
we need to check if J(C^ef) — J(C^'gf), or equivalently, 
J{Cl^) = JiC^^f); this is because J{C^^f) = J{C^^f) due 
to C^^f < C^lf and the definition of C,1f. Due to numerical 



-■ref 



ref 



issues in checking the equality, we introduce a small constant 
Tbost > SO that we treat J{Cl^f) the same as J{C^^f) if 
JiCl.^f)-J{C^^f) < Tcost. Similarly, we treat J{Ci^f) > J{Cf^) 
if J(C|!gf) — J(C|"f) > Tcost- The detailed algorithm is given in 
Algorithm [1] At Step 4, if J(Cjgf) - J(Cf^f) > Tcost, or equiv- 
alendy, J(C;,f) - JiC^^^) > Tcost, we have J(C;,f) > J{C^^). 
Because of the monotonicity property in Proposition |2] we 
know J(Cj,f) > J(C7;,f) > J(C;;i) = JiC^^) for any 
j = i + 1, L. Therefore, the for loop can be terminated, 
and the approximated critical battery capacity is Cjl^j-. 

It can be verified that Algorithm [T] stops after at most L + 1 
steps, or equivalently, after solving at most 



' cap 



1 



optimization problems in Eq. (O. The accuracy of the critical 
battery capacity is controlled by the parameters Tcap, Tcost- 
Fixing St, Tcost, the output is within 



ref 



^cap , C'ref 



'^cap] 



Therefore, by decreasing Tcap, the critical battery capacity can 
be approximated with an arbitrarily prescribed precision. 

Since the function J(Cref) is a non-increasing function of 
Cref, we propose Algorithm |2] based on the idea of bisection 
algorithms. More specifically, we maintain three variables 

Cref < C'ref < ^ref i 



ref 

Input: The optimization problem in Eq. ^ with fixed 
to,T, K, Z,Tc, D under Assumptions [T] and 121 the calculated 
bounds C,'ef, Cref, and parameters 5t, reap, Tcost 

Output: An approximation of CJef 

1: Let Cif = and = C;*; 

2: Solve the optimization problem in Eq. ^ with C^f, and obtain 



Let sign — 1; 
while sign = 1 do 
Let Cref = "•' 2 ; 

Solve the optimization problem in Eq. ^ with Cref, and obtain 



J (Cref 



if J(C^f^ 
Set Cref 
else 

Set Cref = Cref and set J(C, 



< Tcost then 
Cif, and J{Clf) = 



if ai 



C„f < Tcao then 



with J{Cl 



Set sign = 
end if 
end if 
end while 

Output Cref = C'i 



0; 



in which Clf (or Clf) is initialized as Cj^f (or C^^f). We set 
. Due to Proposition |2] we have 



Clf to be 



If J{Cl,) = J{Ci,) (or equivalently, J{Ci,) = J(CS); this 
is examined in Step 7), then we know that C^f < C'^^f < C^^f, 
therefore, we update C^^^ with C^^f but do not updatj] the 
value J{Clf). On the other hand, if J{Clf) > J{Clf), then 
we know that C^^f > C^^f > C^^f, therefore, we update C^j 
with Cj^gf and set J(C/gf) with J{C^^f). In this case, we also 
check if C^^f — C^^f < Tcap: if it is, then output C^^f since we 
know the critical battery capacity is between C^^f and C^^f, 
otherwise, the while loop is repeated. Since every execution 
of the while loop halves the interval [Cj^f , C^.^] starting from 
[Cjgf, Cj"f], the maximum number of executions of the while 

loop is [log2 — '—], and the algorithm requires solving at 
most 

\\og, + 1 

^cap 

optimization problems in Eq. dS). This is in contrast to 
solving \ ^"'^ '~^"' ~\ + 1 optimization problems in Eq. ^ using 
Algorithm Q] 

Remark 4 Note that the critical battery capacity can be 
implemented by connecting batteries with fixed capacity in 
parallel because we only assume that the minimum battery 
charging time is fixed. ■ 

V. Simulations 

In this section, we calculate the critical battery capacity 
using Algorithm |2] and verify the results in Section |lll] via 
simulations. The parameters used in Section HI] are chosen 

'Note that if we update J{C^^f) with J(C^^f), then the difference between 
J(Cj^^j) and J{C"^f) can be amplified when C^^j is updated again later on. 



(a) Starting from July 8, 2011. 



Fig. 3. Time-of-use pricing for the summer season at San Diego, CA 1161 . 



(b) Starting from July 13, 2011. 

Fig. 2. PV output based on GHI measurement at La Jolla, California, where 
tick marks indicate noon local standard time for each day. 



based on typical residential home settings and commercial 
buildings. 

A. Setting 

The GHI data is the measured GHI in July 2010 at La 
Jolla, California. In our simulations, we use rj — 0.15, and 
S = lOm^. Thus Fpv(i) = 1.5 x Gm{t){W). We have two 
choices for to: 

. to is 0000 h local standard time (LST) on Jul 8, 2010, 
and the hourly PV output is given in Fig. |3a) for the 
following four days starting from to (this corresponds to 
the scenario in which there are relatively large variations 
in the PV output); 

. to is 0000 h LST on Jul 13, 2010, and the hourly PV 
output is given in Fig. |2tb) for the following four days 
starting from to (this corresponds to the scenario in which 
there are relatively small variations in the PV output). 

The time-of-use electricity purchase rate Cg{t) is 

. 16.5<f/kWh from 11AM to 6PM (on-peak); 
. 7M/kWh from 6AM to 1 1AM and 6PM to 10PM (semi- 
peak); 

• 6.1(^/kWh for all other hours (off-peak). 

This rate is for the summer season proposed by SDG&E |fT6l . 
and is plotted in Fig. [3] 

Since the battery dynamics and aging are characterized by 
continuous ordinary differential equations, we use St = Ih as 



the sampling interval, and discretize Eqs. (O and (O as 

EB{k + l) = EB{k)+PB{k) , 

[ AC(fc) - ZPb (fc) if Pb (k) < 



AC(fc + l) = 



I AC(fc) 



otherwise 



In simulations, we assume that lead-acid batteries are used; 
therefore, the aging coefficient is Z = 3 x 10^^ |i5J, the unit 
cost for capacity loss is K = 0.15$/Wh based on the cost 
of 150$ /kWh ETI . and the minimum charging time is Tc = 
I2h L22J. 

For the PV DC-to-AC converter and the battery DC-to- 
AC/AC-to-DC converters, we use rjpy — 1]^ — 0.9. It can be 
verified that Assumption |2] holds because the threshold value 
for K iE 

(maxt Cg{t) ~ mint Cg{t))riB _ 



(16.5 - 6.1) X 10-5 X 0.9 



0.3120 



3 X 10-4 

For the load, we consider two typical load profiles: the res- 
idential load profile as given in Fig. SJa), and the commercial 
load profile as given in Fig. SJb). Both profiles resemble the 
corresponding load profiles in Fig. 8 of [15|0 Note that in 
the residential load profile, one load peak appears in the early 
morning, and the other in the late evening; in contrast, in the 
commercial load profile, the two load peaks appear during 
the daytime and occur close to each other For multiple day 
simulations, the load is periodic based on the load profiles in 

Fig. a 

For the parameter D, we use D = 800(M^). Since the 
maximum of the loads in Fig. |4]is around 1000(1^), we will 
illustrate the peak shaving capability of battery storage. It can 
be verified that Assumption [T] holds. 

B. Results 

We first examine the storage size determination problem 
using both Algorithm [T] and Algorithm |2] in the following 
setting (called the basic setting): 

* Suppose the battery is a Li-ion battery with the same aging coefficient as 
a lead-acid battery. Since the unit cost K = 1.333$/Wh based on the cost 
of 1333$ /kWti EU (and K = 0.78$/Wh based on the 10-year projected 
cost of 780$/kWh f21j), the use of such a battery is not as competitive as 
directly purchasing electricity from the grid. 

'However, simulations in 1151 start at 7AM so Fig. |4]is a shifted version 
of the load profile in Fig. 8 in 1151 . 
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P (W): Battery Charging/Discharging Profile 



(a) Residential load averaged at 536.8{iy). 



19 20 21 22 23 24 



Hours of a day 

(b) Commercial load averaged at 485. 6(W). 
Fig. 4. Typical residential and commercial load profiles. 



• the cost minimization duration is T = 24(/i), 
. to is 0000 h LST on Jul 13, 2010, and 

• the load is the residential load as shown in Fig. SJa). 
The lower and upper bounds in Proposition |5] are calculated as 
2667{Wh) and 39269{Wh). When applying Algorithm[I] we 
choose Tcap — lQ{Wh) and Tcost — 10^^, and have L — 3660. 
When running the algorithm, 2327 optimization problems in 
Eq. ([8]) have been solved. The maximum cost Jmax is —0.1921 
while the minimum cost is —0.3222, which is larger than the 
lower bound —2.5483 as calculated based on Proposition |4] 
The critical battery capacity is calculated to be 16089{Wh). 
In contrast, only 



[log 



' cap 



1 = 13 



optimization problems in Eq. (O are solved using Algorithmic 
while obtaining the same critical battery capacity. 

Now we examine the solution to the optimization problem 
in Eq. ^ in the basic setting with the critical battery capacity 
l6089{Wh). We solve the mixed integer programming prob- 
lem using the CPlex solver ll20l . and the objective function is 
J(16089) = -0.3222. Psit), Esit) and C{t) are plotted in 
Fig. 13 a). The plot of C{t) is consistent with the fact that there 
is capacity loss (i.e., the battery ages) only when the battery 
is discharged. The capacity loss is around 2.2{Wh), and 



Ag(to+T) 

Ci-ef 



2.2 

16089 



1.4 X 10" 







which justifies the assumption we make when linearizing the 
nonlinear battery aging model in the Appendix. The dynamic 



o O O O o o o 



E (Wh): Battery Ctiarge Profile 
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(a) 

C ($/Wh): Dynamic Pricing 



Pg(W): Battery Ctiarging/Discharging Profile 
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P JW): Load (solid red curve), and P (W): Net Power Purchase (dotted blue curve) 



-2000 - 
-4000 L- 



(b) 



P (W): Load (solid red curve), and P (W): Net Power Purchase (dotted blue curve) 




(c) 

Fig. 5. Solution to a typical setting in which to is on Jul 13, 2010, T ■ 
24(h), the load is shown in Fig.|4ja), and def = 160S9(Wh). 



pricing signal Cg{t), Psit), fload(i) and Pg{t) are plotted in 
Fig. HJb). From the second plot in Fig. HJb), it can be observed 
that the battery is charged when the time-of-use pricing is low 
in the early morning, and is discharged when the time-of-use 
pricing is high. From the third plot in Fig. IHb), it can be 
verified that, to minimize the cost, electricity is purchased 
from the grid when the time-of-use pricing is low, and is 
sold back to the grid when the time-of-use pricing is high; 
in addition, the peak demand in the late evening (that exceeds 
D = 800(VF)) is shaved via battery discharging, as shown in 
detail in Fig. EJc). 

These observations also hold for T = 48(/i). In this case, we 
change T to be 48(/i) in the basic setting, and solve the battery 
sizing problem. The critical battery capacity is calculated to be 
Qrf ~ 16096(VF/i). Now we examine the solution to the opti- 
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P (W): Battery Charging/Discharging Profile 
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C ($/Wh): Dynamic Pricing 
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(b) 







Fig. 6. Solution to a typical setting in which to is on Jul 13, 2010, T ■ 
48(/i), the load is shown in Fig.EJa), and Cef = 16096(Wh). 



mization problem in Eq. (HJ with the critical battery capacity 
16096(VF/i), and obtain J(16096) = -0.6032. Psit), Esit), 
C{t) are plotted in Fig. ISJa), and the dynamic pricing signal 
Cg{t), Psit), Pioad(i) and Pg{t) are plotted in Fig.l^b). Note 
that the battery is gradually charged in the first half of each 
day, and then gradually discharged in the second half to be 
empty at the end of each day, as shown in the second plot of 
Fig.lSa). 

To illustrate the peak shaving capability of the dynamic 
pricing signal as discussed in Remark [T] we change the load 
to the commercial load as shown in Fig. Ub) in the basic 
setting, and solve the battery sizing problem. The critical 
battery capacity is calculated to be C^^f — l3352{Wh), and 
J(13352) — —0.1785. The dynamic pricing signal Cg{t), 
Psit), Pioad(i) and Pg{t) are plotted in Fig. |7] For the 
commercial load, the duration of the peak loads coincides 
with that of the high price. To minimize the total cost, during 
peak times the battery is discharged, and the surplus electricity 
from PV after supplying the peak loads is sold back to the 
grid resulting in a negative net power purchase from the grid, 
as shown in the third plot in Fig. [T] Therefore, unlike the 
residential case, the high price indirectly forces the shaving of 
the peak loads. 

Now we consider settings in which the load could be either 
residential loads or commercial loads, the starting time could 
be on either Jul 8 or Jul 13, 2010, and the cost optimization 
duration can be 24(ft,), 48(/i), 96{h). The results are shown in 
Tables U and M In Table H] is on Jul 8, 2010, while in 



Fig. 7. Solution to a typical setting in which to is on Jul 13, 2010, T : 
24(h), the load is shown in Fig.|4jb), and Crf = 13352(Wh). 



Table HH to is on Jul 13, 2010. In the pair (24, i?), 24 refers 
to the cost optimization duration, and R stands for residential 
loads; in the pair (24, C), C stands for commercial loads. 

We first focus on the effects of load types. From Tables J] 
and ini the commercial load tends to result in a higher cost 
(even though the average cost of the commercial load is 
smaller than that of the residential load) because the peaks of 
the commercial load coincide with the high price. When there 
is a relatively large variation of PV generation, commercial 
loads tend to result in larger optimum battery capacity C^^f as 
shown in Table |I] presumably because the peak load occurs 
during the peak pricing period and reductions in PV production 
have to be balanced by additional battery capacity; when there 
is a relatively small variation of PV generation, residential 
loads tend to result in larger optimum battery capacity C^^f as 
shown in Table HIl If tg is on Jul 8, 2010, the PV generation 
is relatively lower than the scenario in which to is on Jul 13, 
2010, and as a result, the battery capacity is smaller and the 
cost is higher. This is because it is more profitable to store 
PV generated electricity than grid purchased electricity. From 
Tables |T] and |II] it can be observed that the cost optimization 
duration has relatively larger impact on the battery capacity 
for residential loads, and relatively less impact for commercial 
loads. 

In Tables U and |ll] the row Jmax corresponds to the cost in 
the scenario without batteries, the row J{C^^f) corresponds to 
the cost in the scenario with batteries of capacity C^^f, and the 
row 5'ay/wgi'l corresponds to J^ax — •^(C'ref)- For Table U we 
can also calculate the relative percentage of savings using the 



formula ' j ^ , and get the row Percentage. One obser- 
vation is that the relative savings by using batteries increase as 
the cost optimization duration increases. For example, when 
T — 96(h) and the load type is residential, 16.10% cost can be 
saved when a battery of capacity 18320{Wh) is used; when 
T = 96{h) and the load type is commercial, 25.79% cost 
can be saved when a battery of capacity 15747{Wh) is used. 
This clearly shows the benefits of utilizing batteries in grid- 
connected PV systems. In Table [III only the absolute savings 
are shown since negative costs are involved. 

'"Note that in Table ITU part of the costs are negative. Therefore, the word 
"Earnings" might be more appropriate than "Savings". 



12 



TABLE I 

Simulation Results for to on Jul 8, 2011 





(24, R) 


(48, R) 


(96, R) 


(24, C) 


(48, C) 


(96, C) 




6891 


7854 


18320 


8747 


13089 


15747 




0.8390 


1.5287 


1.9016 


0.9802 


1.7572 


2.3060 


J max 


0.9212 


1.7027 


2.2666 


1.1314 


2.1232 


3.1076 


Savings 


0.0822 


0.1740 


0.3650 


0.1512 


0.3660 


0.8016 


Percentage 


8.92% 


10.22% 


16.10% 


13.36% 


17.24% 


25.79% 


TABLE II 

Simulation Results for to on Jul 13, 201 1 




(24, R) 


(48, R) 


(96, R) 


(24, C) 


(48, C) 


(96, C) 




16089 


16096 


70209 


13352 


17363 


32910 




-0.3222 


-0.6032 


-0.9930 


-0.1785 


-0.3721 


-0.5866 




-0.1921 


-0.3395 


-0.4649 


0.0181 


0.0810 


0.3761 


Savings 


0.1301 


0.2637 


0.5281 


0.1966 


0.4531 


0.9627 



VI. Conclusions 

In this paper, we studied the problem of determining the 
size of battery storage for grid-connected PV systems. We 
proposed lower and upper bounds on the storage size, and 
introduced an efficient algorithm for calculating the storage 
size. In our analysis, the conversion efficiency of the PV 
DC-to-AC converter and the battery DC-to-AC and AC-to- 
DC converters is assumed to be a constant. We acknowledge 
that this is not the case in the current practice, in which 
the efficiency of converters depends on the input power in 
a nonlinear fashion |5 |. This will be part of our future work. 
In addition, we would like to extend the results to distributed 
renewable energy storage systems, and generalize our setting 
by taking into account stochastic PV generation. 



Appendix 

Derivation of the Simplified Battery Aging Model 

Eqs. (11) and (12) in |5| are used to model the battery 
capacity loss, and are combined and rewritten below using the 
notation in this work: 



C(t + St)-C{t) ^ -C,-efXZx( 



Esit) EB{t + 5t) 
C{t) C{t + St) 



) . (15) 



If 5t is very small, then C{t + St) « C{t). Therefore, 



EB(t+St) 
C{t+St) 



-PB{t)St 

C(t) 



C{t) 

If we plug in the approximation, divide 



5t on both sides of Eq. dTsl l. and let 5t goes to 0, then we 
have 



dC{t) 
dt 



C[t) 



(16) 



This holds only if Psit) < as in IS), i.e., there could be 
capacity loss only when discharging the battery. 

Since Eq. (fT6] l is a nonlinear equation, it is difficult to solve. 
Let AC{t) = (7ref — C{t), then Eq. ( fTSI l can be rewritten as 



dACjt) 
Jt 



-Z 



PB{t) 

_ AC{t) 



(17) 



If t is much shorter than the life time of the battery, then the 



percentage of the battery capacity loss 



AC(t) 
C,,, 



0. Therefore, Eq. (fTTI i can be simplified to the following linear 
ODE 

when Psit) < 0. If Psit) > 0, there is no capacity loss, i.e., 

dC{t) 



dt 







is very close to 
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